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ABSTRACT Systemic lupus erythematosus (SLE) is an 
autoimmune multisystem inflammatory disease characterized 
by the production of pathogenic autoantibodies. Previous 
genetic studies have suggested associations with HLA Class II 
alleles, complement gene deficiencies, and Fc receptor poly- 
morphisms; however, it is likely that other genes contribute to 
SLE susceptibility and pathogenesis. Here, we report the 
results of a genome-wide microsatellite marker screen in 105 
SLE sib-pair families. By using multipoint nonparametric 
methods, the strongest evidence for linkage was found near the 
HLA locus (6pll-p21) [D6S257, logarithm of odds (lod) = 
3.90, P = 0.000011] and at three additional regions: 16ql3 
(D16S415, lod = 3.64, P = 0.000022), 14q21-23 (D14S276, 
lod = 2.81, P = 0.00016), and 20pl2 (D20S186, lod = 2.62, P = 
0.00025). Another nine regions (lp36, lpl3, lq42, 2pl5, 
2q21-33, 3cent-qll, 4q28, llplS, and 15q26) were identified 
with lod scores ^1.00. These data support the hypothesis that 
multiple genes, including one in the HLA region, influence 
susceptibility to human SLE. 



Systemic lupus erythematosus (SLE) is an idiopathic autoim- 
mune disease in which self-reactive autoantibodies cause dis- 
ease either by directly binding to self-antigens (e.g., antiphos- 
pholipid syndrome, immune cytopenias) or following the 
deposition of antibody-antigen immune complexes in blood 
vessels leading to vasculitis and tissue damage (1). The esti- 
mated prevalence of SLE in the U.S. is ^45/100,000, with the 
peak incidence in women ages 20-40 (2). The familial clus- 
tering of SLE, together with higher rates of concordance in 
monozygotic compared with dizygotic twins, suggests an im- 
portant role for genetic predisposition to systemic lupus (3-6). 

Although the clinical manifestations of systemic lupus are 
heterogeneous, SLE can be distinguished from a number of 
related syndromes (drug-induced lupus, discoid lupus, sub- 
acute cutaneous lupus) and other diseases by careful clinical 
examination and laboratory testing. Diagnostic criteria for 
SLE have been developed by the American College of Rheu- 
matology that are both highly specific and sensitive (7, 8). 
Because of the likelihood that multiple genes are responsible 
for susceptibility to SLE, and with the sibling relative risk ratio 
(A s ) (9) for SLE ~ 10-20 (10), we elected to use a sibling 
(sib)-pair family approach to map potential genetic loci. 

MATERIALS AND METHODS 

Families. The recruitment of families for this study has 
been described (11). All patients met the 1997 American 
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College of Rheumatology revised criteria for the diagnosis of 
SLE (7, 8). 

Samples and Genotyping. Genomic DNA was isolated from 
peripheral blood cells by using standard conditions. Genotyp- 
ing was performed by using an Applied Biosystems fiuores- 
cently labeled human linkage mapping set (version 1.0, panels 
1-18; version 2.0, panels 19-27). PCR (32 cycles) was per- 
formed on an ABI 877 Catalyst robotic workstation [5 yl 
reactions; 8 ng of genomic DNA, 2.5 mM MgCb, 0.2 mM 
dNTPs (Pharmacia), 1.65 pmol of 5' and 3' primers, 0.2 units 
of Amplitaq Gold DNA polymerase (Perkin-Elmer) in IX 
PCR Buffer II (Perkin-Elmer)]. 

Pooled amplification products were electrophoresed 
through 5% polyacrylamide gels (FMC) for 2 hr at 3,000 V by 
using an ABI 377 DNA Sequencer. Semiautomated fragment 
sizing was performed by using genescan 2.1 software (ABI) 
followed by allele calling with genotyper 2.0 software (ABI). 
Each genotype was reviewed manually by two members of the 
research team to confirm the accuracy of allele calling. Poor 
performing markers (a total of 17) were excluded from the 
analysis. The overall data drop-out rate for the ** 127,500 
genotypes analyzed was <1.5%. 

Data Analysis. Nonparametric multipoint analysis was per- 
formed with genehunter plus (12), a modified version of 
GENEHUNTER (13), by using the "all" statistic. Allele frequen- 
cies for the parameter file were generated from the unaffected 
parental genotypes for the cohort. Marker map positions were 
obtained from the sex-averaged maps compiled by J. Weber 
(Marshfield Clinic, Marshfield, WI) (www.marshmed.org/ 
genetics/). 

RESULTS AND DISCUSSION 

Families were recruited by advertising for "Sisters with Lupus" 
(11), and, as a result, our sample was highly enriched for female 
patients (female:male 219:1), compared with the population 
estimate for lupus of ^90% female. Selected clinical and 
demographic data of the SLE patients in this study are shown 
in Table 1. 

The sample studied in this report comprises the first 105 SLE 
sib-pair families collected (102 families with 2 SLE sibs, 2 
families with 3 SLE sibs, and 1 family with 4 SLE sibs). All 
available parents (including three affected) were collected, 
and, in the absence of parents, an unaffected sibling was 
sampled to assist in reconstruction of parental genotypes. In 
three families an additional affected-first degree relative was 
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Table 1. Clinical and demographic features of 220 SLE patients 



Aoe at diapnosis **~SD 


31 ■*■ 11 vears 


Duration nf Hi^pn^p 




Sex female'male 


219*1 


Fthnicitv % 






80 


T-Ticnnnif 


8 


A f rifa n - A m p rica n 


5 


Asian 


3 


N^iypH hpritncp 


4 


T ahoratorv clinical features* % 




r\l > n |JU5 1 live 


98 


Anti-d^RNA-nositive 


46 


Arthritis 


85 


OM11 IIIVUIVCIIICJII 


Q1 


i icur ilia 




Hematologic 


47 


Renal disease 


30 


CNS lupus 


22 


Pericarditis 


18 


More than one miscarriage 


13 


Medication history, % 




Corticosteroids 


77 


Antimalarials 


65 


Cytotoxic drugs 


28 


Intravenous steroids 


19 



*Data represents the percentage of SLE patients having the indicated 
laboratory/clinical features and medication history at any time 
during the course of their disease. 

1"Four individuals negative for ANA tested positive for anti-dsDNA 
antibodies and otherwise fulfilled criteria. These data are comparable 
with those described in other large series of SLE patients (31). ANA, 
antinuclear antibodies; dsDNA, double-stranded DNA; CNS, central 
nervous system. 

included in the analysis. The final study cohort totaled 220 SLE 
patients and 155 unaffected parents or sibs. The racial com- 
position of the 105 families studied was as follows: 84 Cauca- 
sian, 8 Hispanic, 6 African-American, 3 Asian, and 4 of mixed 
heritage. 

Family members were genotyped with 341 highly polymor- 
phic markers across the 22 autosomes, at an average inter- 
marker distance of 9.7 centimorgans (cM). Mendelian inher- 
itance was confirmed in all families, and nonparametric mul- 
tipoint linkage analysis was performed on the collected marker 
data by using genehunter (13), with extensions to include 
calculation of appropriate logarithm of odds (lod) scores and 
support intervals (12). 

Results of the multipoint analysis are shown in Fig. 1, with 
Zi r , the test statistic generated by genehunter plus, and the 
calculated lod scores for each marker plotted for each chro- 
mosome. We considered any region with a Z\ T > 2.0 as 
"potentially interesting," and applied the criteria of Lander 
and Kruglyak (14) to further define regions of significant (Z >: 
4.1, lod > 3.6, P < 0.00002) or suggestive (Z > 3.2, lod > 2.2, 
P < 0.0007) linkage. Overall, 25 of the 341 markers tested (7%) 
gave Zi r scores >2.1 and lod scores >1.00. Strikingly, 16 of the 
25 positive markers clustered into 4 distinct genomic intervals 
(Fig. 1 and Table 2). Two of these intervals contained a marker 
that met criteria for significant linkage: 6pll-p21, mapping 
near the HLA region (D6S257, Z, r = 4.24, lod = 3.90, P = 
0.000011), and 16ql3 (D16S415, Zi r = 4.09, lod = 3.64, P = 
0.000022). Two other regions fulfilled criteria for suggestive 
linkage: 14q21-q23 (D14S276, Z ir = 3.60, lod = 2.81, P = 
0.00016) and 20pl2 (D20S186, Z, r = 3.48, lod = 2.62, P = 
0.00025). Nine additional chromosomal regions were identi- 
fied by single markers with Z\ T scores >2.1 and accompanying 
lod scores 5:1.00 (Table 2). In all of the most suggestive areas, 
these multipoint results were confirmed by single point anal- 
yses (data not shown). 
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Because of sample size considerations, we stratified the 
data by ethnic group only for the 84 Caucasian families in the 
overall sample of 105 families. The results of this analysis are 
shown in Table 3. Lod scores dropped in 10 of the top 13 
identified potential susceptibility intervals when non- 
Caucasian families were eliminated from the analysis, sug- 
gesting that families of all ethnic groups contributed to the 
evidence for genetic linkage in these regions. In contrast, 
three intervals (4q28, llpl5, and 15q26) were characterized 
by an improvement in lod scores when only the Caucasian 
families were considered. 

The strongest evidence for linkage was found at 6pll-p21, 
with the best markers (D6S426 and D6S257) mapping just 
centromeric to the HLA region (located at 6p21.3). There is 
a long history of interest in the role of the major histocom- 
patibility complex (MHC) in many autoimmune diseases, 
including IDDM, familial psoriasis, multiple sclerosis, rheu- 
matoid arthritis, and SLE. Genetic associations have been 
shown for HLA-DR2 or HLA-DR3 alleles in SLE patients 
from several ethnic groups (relative risks ranging from 2.0 to 
3.4) (15). However, in large lupus-prone families, HLA 
alleles do not necessarily segregate with disease (16). A 
stronger association of HLA alleles (DQ in particular) has 
been demonstrated for specific autoantibody subsets in SLE 
(17). Non-HLA genes within the MHC, notably the tumor 
necrosis factor a gene and the complement components C2 
and C4, also have been implicated in SLE (17). Further work 
will be necessary to determine whether this 6p locus repre- 
sents a polymorphic HLA locus or a linked gene within or 
near HLA, 

A recent comparison of the linkage results from 23 
independent genome screens in various human and experi- 
mental animal autoimmune or inflammatory diseases iden- 
tified 18 "clusters" of non-MHC candidate human autoim- 
mune loci (18). Of interest, two of the non-MHC intervals 
identified in this screen of SLE families (16ql3 and llpl5) 
are located within clusters implicated in other human auto- 
immune diseases. The 16q cluster is a large interval of «35 
cM that includes potential susceptibility loci for Crohn's 
disease (19), Blau syndrome (20), psoriasis (21), IDDM (22), 
and asthma (23). Cluster lip is a narrow interval (^2-3 cM) 
identified in asthma (23), multiple sclerosis (24), and IDDM 
(22). SLE is not known to associate strongly with any of the 
diseases that define these clusters, but these data suggest the 
possibility that there is a sharing of genetic predisposition to 
multiple autoimmune diseases or, alternatively, that each of 
these intervals may contain closely linked autoimmune 
susceptibility genes. 

The results of an independent genome scan of SLE families 
performed by Moser et al (25) at the Oklahoma Medical 
Research Foundation support many of the findings of this 
screen. Of the top 13 intervals identified in this sib-pair family 
screen, 4 (20pl2, lpl3, lq42, and 4q28 with lod scores ranging 
from 2.62 to 1,46) also were identified in the Oklahoma 
Medical Research Foundation screen of mostly larger families 
multiplex for SLE. The 20pl2 region shows suggestive evi- 
dence for linkage in this study (D20S186, Z )r = 3.48, P = 
0.00025, lod = 2.62) and was also one of the strongest loci 
identified in the Oklahoma screen (25). This interval maps 
centromeric to a recently identified susceptibility interval for 
psoriasis (21). The locus in the lq42 region, originally identi- 
fied by Tsao et al in a candidate approach in 52 SLE sib-pairs 
(26), now has been identified in three independent SLE 
screens [refs. 25 and 26 and this study (D1S235, Zi r = 2.64, P = 
0.0041, lod = 1.51)]. Despite these similarities, several of the 
strong intervals identified in this screen (6p, 14q and 16q) did 
not show evidence for linkage in the Oklahoma Medical 
Research Foundation study. Also, in our primarily Caucasian 
sample, we find little evidence for linkage at lq22-23. This 
region showed evidence for significant linkage in the Okla- 
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Fig. 1. Nonparametric multipoint linkage analysis in 105 SLE sib-pair families. Shown are the analyses of chromosomes 1-22, plotting both Z\ T 
(solid line) and lod (dotted line) scores. Subjects were genotyped with 341 polymorphic markers at an average interval of 9.7 cM (longest, 26.4 
cM; shortest, 0.5 cM). The average value of the information statistic, which is a measure of the percent inheritance information extracted at each 
marker (13), was 0.78, indicating that, overall, the markers were highly informative. 



homa Medical Research Foundation study and harbors the in African-American SLE patients (27). These differences in 
various Fc receptors, which are associated with renal disease the genome screen results may reflect "false positive" evidence 
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Table 2. Summary of potential SLE susceptibility loci 



Zir range 



Interval* 


z, r 


P value 


lod score* 


Positive markers 


Map position, cM* 


IC§ 


6pll-p21 


3.69 


0.00011 


2.96 


D6S426 


60.44 


0 77 




4.24 


0.000011 


3.90 


D6S257 


79.92 


0 79 




2.19 


0.014 


1.04 


D6S462 


99.01 


0 63 


16ql3 


3.81 


0.000069 


3.15 


D16S3136 


62.11 


0.77 




a no 


u.uuuuzz 


J. 04 


D10o41j 


07. oz 


0.73 




2.42 


0.0077 


1.27 


D16S503 


83.55 


0 69 


14q21-23 


2.94 


0.0016 


1.87 


D14S288 


47.51 


0.70 




3.60 


0.00016 


2.81 


D14S276 


56.36 


0.77 




9 8ft 
Z.oU 


U.UUZO 


1 71 
1. /I 




aa 1 o 


A OI 




2.99 


0.0014 


1.94 


D14S258 


76.28 


0.78 




2.54 


0.0055 


1.40 


D14S74 


87.36 


0.50 




9 90 

Z.Z7 


ft m 1 

U.Ul I 


1 14 


r»14Q9ftft 




U.0/ 




% 48 


n nnn9s 

U.UUUZJ 


9 69 
z.oz 


nOflA 1 R6 


19 1ft 


U.oU 




3.21 


0.00066 


2.23 


D20S112 


39.25 


0.80 




2.38 


0.0087 


1.23 


D20S107 


55.74 


0.87 




9 77 
Z. / / 


n nn9ft 


1 fn 

i.O/ 


ronci 1 o 


£1 '7'7 

Oi. / / 


C\ TO 


zpi j 


9 7ft 

Z. /o 


ft ftft97 


l.Oo 




ou.oy 


f\ OC 

U.oj 


1r»1 1 


Z.OO 


ft ftftlQ 


1 ^i 


F*1 C9^9 


1 T7 


n ta 
U.7o 


1n49 


9 64 


ft ftft41 


i ^1 


ni coic 
U 1 ozj j 


ZD4.04 


U.67 


4q28 


2.60 


0.0047 


1.46 


D4S424 


144.46 


0.74 


3cent-qll 


2.39 


0.0084 


1.24 


D3S1271 


117.76 


0.57 


llpl5 


2.34 


0.0096 


1.19 


D11S922 


2.11 


0.86 


2q21-33 


2.27 


0.012 


1.12 


D2S151 


152.04 


0.76 


15q26 


2.22 


0.013 


1.07 


D15S127 


86.81 


0.86 


lp36 


2.15 


0.016 


1.00 


D1S234 


55.10 


0.74 



>4.0 



3.0-4.0 



2.0-3.0 



*Chromosome locations were determined from the maps available at www.gdb.org/and http://cedar.genetics.soton.ac.uk/pub/. 
tlod - Z, r 2 /2lnl0 (12). 

*Map positions obtained from the sex-averaged maps compiled by J. Weber (www.marshmed.org/genetics/). 
^Information content (13) for the marker designated in each interval. 



for linkage at some loci (14) or may result from ethnic and/or 
genetic heterogeneity in the samples. 

In summary, using an affected sib-pair family approach, we 
have identified four potential SLE loci with lod scores >2.6 
and an additional nine loci with lod scores ranging from 1.00 
to 1.68. The four most interesting regions (6pll-p21-HLA, 
14q21-23, 16ql3, and 20pl2) show surprisingly strong evi- 
dence for linkage given the sample size of 114 affected sib pairs 
and a total of 127 affected relative pairs. Two of the non-MHC 

Table 3. Comparison of lod scores obtained for Caucasian SLE 
sib-pair families versus entire family cohort 



Loci 


Caucasian families* 


All families t 


Difference 




lod score* 




6pll-p21 


2.50 


3.90 


-1.40 


16ql3 


3.36 


3.64 


-0.28 


14q21-23 


2.56 


2.81 


-0.25 


20pl2 


2.09 


2.62 


-0.53 


2pl5 


1.50 


1.68 


-0.18 


lpl3 


1.27 


1.53 


-0.26 


lq42 


0.46 


1.51 


-1.05 


4q28 


2.00 


1.46 


+0.54 


3cent-qll 


0.87 


1.24 


-0.37 


llpl5 


1.41 


1.19 


+0.22 


2q21-33 


0.61 


1.12 


-0.51 


15q26 


2.09 


1.07 


+ 1.02 


lp36 


0.41 


1.00 


-0.59 



intervals identified in this screen (16ql3 and llpl5) have been 
implicated in other human autoimmune diseases. These data 
support results from genetic studies in lupus-prone mouse 
strains suggesting that multiple genes are responsible for 
conferring susceptibility to SLE (28-30). It will now be im- 
portant to confirm these findings in additional cohorts of 
sib-pair families and to initiate efforts to further narrow these 
intervals in preparation for gene identification. 

We are grateful to all of the patients and their family members who 
participated in this study and the many physicians who referred 
families and verified diagnoses. Thanks to J. Harley and K. Moser for 
sharing data before publication. Supported by grants from the Min- 
nesota chapters of the Arthritis Foundation and the Lupus Foundation 
of America, the Texas Gulf Coast Chapter of the Lupus Foundation 
of America, and National Institutes of Health Grants R01 AR/ 
AI43271, SLE SCOR P50 AR45231, and KO8-AR01986. P.M.G. is 
supported by a Physician-Scientist Award from the University of 
Minnesota Department of Medicine. 



*84 Caucasian families. 
tl05 total sib -pair families, 
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